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The problem of propagation of plane electromagnetic waves polarized in the plane of incidence in 
a plane-stratified isotropic plasma (regarded as a dielectric) is dealt with by reducing it to the solution 
of an ordinary differential equation which is singular at the level of plasma resonance. The absorption 
of power in the resonant layer is calculated approximately in two cases involving a linear profile. In 
the first the relative dielectric constant varies linearly in a thin transition layer between two homoge- 
neous regions. In the second, the dielectric constant is linearly varying over a large interval. 

1. Introduction 

The problem of the reflection of a plane electromagnetic wave from a plane stratified plasma 
presents considerably more difficulty in the "vertical" polarization case (i.e., in the case of electric 
field in the plane of incidence) than in the "horizontally" polarized case [Budden, 1961]. In the 
latter, the ordinary differential equation which one obtains after eliminating the variables trans- 
verse to the stratification gradient exhibits a turning point at some level of the stratification, while 
in the former case, in addition to the turning point, the differential equation has a singularity. This 
singularity occurs at the stratum of plasma resonance. 

We will assume that the plasma is an isotropic dielectric medium, with the dielectric constant 
related to the electron density N(z) by the well-known formula 

e = e (l — Ne 2 /m€ (D 2 ), (i) 

where e and m are the charge and mass of the electron, and with permeability equal to that of free 
space. We assume further that a time harmonic wave (with time dependence e> w ) propagates in 
the medium. The physical inconsistencies arising from power loss at the level where e = in a 
lossless medium are resolved by appealing to the infinite duration of the transient process [Ginz- 
burg, 1961], to nonlinearities which must occur when the fields become large [Budden, 1961], 
to coupling to plasma waves [Denisov, 1957], or to the impropriety of the dielectric model for the 
plasma, or they are avoided by considering a lossy plasma [Wait and Walters, 1963]. We will 
not concern ourselves here with these explanations. We shall merely obtain an approximate solu- 
tion of the differential equation for a linearly varying electron density. The results will then be 
compared K) those obtained previously by numerical methods [Budden, 1961]. 

1 Taken from the dissertation submitted by Peter Hirsch to the Faculty of the Brooklyn Polytechnic Institute in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy, 1964. 
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FIGURE 1. Physical configuration. 



2. Linear Profile 

In a plane stratified dielectric medium (see fig. 1), a plane wave incident in the xz-plane, 
polarized in the plane of incidence constitutes an E-mode with H y , E Xy and E z as the only non- 
vanishing field components. After separating the field dependence on the x coordinate through 
the factor exp(— jk x x), we obtain the ordinary differential equation satisfied by H y : 



d 2 H 

dz 2 ' 



dejdz dH 

e dz 



+ ((o 2 fjLe-k 2 .)H = 



(2) 



where the subscript y was omitted from the magnetic field H since the other components are zero, 
and e and /jl denote the dielectric constant and permeability of the stratified medium. We see that 
the differential equation is singular at the zeros of e. 

We shall consider in this paper two special cases, both characterized by a linearly varying 
dielectric constant. First we shall consider a transition from a region of constant negative di- 
electric constant to free space, assuming that the transition layer is thin. Second, we shall ex- 
amine the relationship between the asymptotic forms of the solution of the problem of an infinite, 
linearly stratified medium on the two sides of the zero of the dielectric constant. Although the 
latter problem has no physical significance per se, it yields the connection formula necessary for 
the asymptotic solution of the problem of propagation in a medium whose stratification is both more 
complicated and more realistic. 

If we introduce a linearly varying dielectric constant 



e(z) = 



' 8 ' 



(3) 



where e is the dielectric constant of free space, k the free space wave number and 8 is the slope 
parameter, and then transform the independent variable 



x = 8 1/3 A: z, 



the differential equation (2) takes the form 






(4) 



(5) 



B = 8^(k x lk f. 



The differential equation is singular at the origin. When B = 0, this singularity is only superficial. 
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In this case the direction of incidence is parallel to the stratification gradient, and the solutions of 
the differential equation are derivatives of the Airy functions [Budden, 1961]. We can express 
the general solution of (5) in terms of power series [Forsterling, 1949] 



H = y[y*;i(x) -f v>{x)\ 



(6) 



whe 



Vl(x) = 2 a nX n * 



(7) 



a =l, «i=0, 



B 



B 



a n = ~ 



q B -3 + fia B -2 



v 2(x)=-v 1 (x)lnx + y d p n x n 



(8) 



j3 =l, j8i=j8 2 = 0, 



fin 



= flj8 M -2+j8 w - 3 -fi(yi-l)a w - 
n(n — 2) 



If we choose to normalize our solution in such a way that //(0) = 1, as is convenient in this 
problem, we simply set y= 1 and retain y as the arbitrary constant to be determined from the 
boundary conditions of the problem. 

To make use of this solution in a physical problem we must resolve the question of the proper 
sheet of the logarithmic function; this can be done by considering the medium to have small losses, 
or by requiring the singularity to be a power sink rather than a power source. 



3. The Thin, Rapidly Varying Layer 

If the linearly stratified layer occurs as a transition between two homogeneous regions (cf. 
fig. 2) of relative dielectric constants 1 and — 63, we can solve the reflection problem by making 
use of the power series (7) and (8) directly. In order that the series converge rapidly, we shall 
assume that the layer is thin; i.e., that the values of the dimensionless variable x at both boundaries 
of the layer be small compared to unity in magnitude. 

Recognizing that in region III the wave is exponentially decaying, we apply, at #3 = 8 2/3 €3 the 
boundary condition 



dH/dx = - (V€ 3 8 2 ' 3 + fl) //. 



(9) 



«/«o 



FIGURE 2. The thin, rapidly varying layer. 
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FIGURE 3. Absorption in a thin transition layer. 



ANGLE OF INCIDENCE 



Applying this boundary condition permits us to evaluate y in the normalized general solution (6) 

(10) 



Ve^ + H 



2e 3 8" 3 A;o 



^8 + 0(8*)] 



2k 



where only the linear term in 8 was retained in the bracket since the present discussion applies 
to the thin, rapid transition only. By matching boundary conditions at X\, we obtain the reflection 
coefficient in region I, T, again to first order in 8 



r _ - eg Vkfl^+j ygj +g r , 2e 3 Vjkl - *» Xc 3 &g + k\) (Q +jP) 



e 3 y/ki-ki+jVe,ki+ki 



€*(**-*») + **» + *» 



+ 0(8 2 



Ves + AJ/A? 



2 O 6 3 Ve3 + H/l 

The power absorption coefficient is the quantity of physical interest here; we obtain 



(id 

(Ha) 
(lib) 



_4 7r8(4/^)6 2 Vl -**/*£ 
e!(l-^Mg) +e 3 + &!/*! 



(12) 



One can obtain this result most conveniently by considering the energy transport associated with 
the wave propagated in the inhomogeneous medium. The dependence of this fractional power loss 
(divided by the factor 47r8e 3 ), on angle of incidence is shown on figure 3. Since 8e 3 remains con- 
stant if the transition layer thickness remains constant, these curves represent the dependence 
of loss on angle of incidence and plasma density for a fixed transition layer thickness. We see 
that maximum absorption occurs at an angle of incidence of about 60° and €3 ~ 1.7. 



4. Connection Across the Singular Point 

To obtain connection formulas for asymptotic forms of the solution of (5) we resort to Langer's 
method [Langer, 1949] . In order to eliminate the first derivative term, we transform the depend- 
ent variable 



that 



U(x) = x~ x m(x) 

d 2 U/dx 2 - [x + B + 3/(4x 2 )] U = 0. 
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(13) 



(14) 



The comparison equation must have identical singular and constant terms and an appropriate 
monotonic term in the coefficient of U. The equation which will be used here has a quadratic 
rather than linear term. This has the consequence that, in order to accomplish the change of 
sign of that term at the origin, we must employ different comparison equations for x > and x < 0. 
The procedure to be followed is thus going to be slightly more complicated than the usual phase 
integral method. In the first step, we derive the approximate behavior of the solution at x = 

from the boundary condition at x > + °°. Normally we would, say, evaluate the logarithmic 

derivative (or impedance) at x = 0; in the present case this is zero for all solutions of the differen- 
tial equation. As was mentioned above the information needed in our case is the value of the 
coefficient y. Once the value of y is obtained, we make use of it as a boundary condition at the 

origin to obtain the asymptotic solution for x > — °°. 

In the region < x < °°, we will obtain an approximation to U in terms of the solution of the 
equation 



d 2 V+ld£ 2 - [a( 2 + B + 3/(4£ 2 )F+ = . 

The relationship between f and x is the usual one in phase integral procedures: 

J vV+£ + 3/4(*') 2 dx' = fVa(£') 2 + £ + 3/[4(f) 2 ] d? 
Jo Jo 



(15) 



(16) 



which, for small x, reduces to 



k*)«* [i +!**+()(*<)]• 



(16a) 



The solution of (15) which decays exponentially as £ > °° is related to the Whittaker function 

[Buchholz, 1953]r x , M by 



r^c^r^ifvaf), 



4Va 2 



(17) 



One can verify readily that (16) and (17) define a function which has the same asymptotic character 
for large positive x as U, has the same behavior at the origin, and satisfies a differential equation 
which differs in the coefficient of the term in U by a function which is everywhere small compared to 
unity [Langer, 1949]. The parameter a has been left undetermined, so far. Let us compare the 
series expansions for the exact solution of (5) and the approximation in terms of Whittaker func- 
tions. We have 



H(x) = H(0) 



and 



h + (x) = 



c, 



l+-x 2 ln x + yx 2 + 0(x 3 ) 



B 



r(l+B/4 Va) 



1 +- x 2 In x + y + {B)x 2 + 0(* 3 ) 



(18) 



(19) 



fhere 



y+(B)= 



n 



<//(lfi/4 V^)+2C-| + | In a 



'2 ' 



(20) 



T is the Gamma function, i// its logarithmic derivative, and C ( > Eulers constant. 
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Figure 4. Three approximations to y(B). 



Figure 5. Plot of 0(B). 



The constant y in the exact solution is a function of B. If we regard a in (20) as a constant 
and evaluate it so that y+(0) is exact then it turns out that y+{B) is a good approximation to y(B) 
for all values of B. We obtain y+(0) (in terms of Airy functions) from the known solution of the 
normal incidence problem [Budden, 1961]. Figure 4 shows a comparison between the function 
y+{B) thus obtained and another approximation in terms of Hankel functions used by Denisov 
[1957] . The zero order approximation is y ; the first order approximation, obtained by converting 
the differential equation into an integral equation and iterating once, is yi. We see that for large 
values of B, y , y u and y+ are nearly identical. The agreement between yi and y + is good for all 
B except at B = 0, where yi is known to be in error and y+ is exact. 

Having calculated y, we can now determine the asymptotic form of the solution for jc — > — oo. 
In this case the comparison equation is 



dW-ldr] 2 - 



-W + B + —)* 



The transformation of variables is again obtained by equating the phase integrals 

( y V-y' +£ + 3/ [4(y') 2 ]dy'= f V-/3(t0 2 + B + 3/[4(V) 2 ] dq' 
Jo Jo 



(21) 



(22) 



fhere 



y=\x\. 



Here the constant /3 is determined, and not arbitrary as in the preceding case, by the requirement 
that (22) must map the turning point of (14) into the turning point of (21). Thus, if y and 170 are 
the zeros of the radicands in (22), /3 is determined by 



fyo fi?o 

V-y + £ + 3/(4y 2 )dy= V-^ 2 -h£+3/(4r ? 2 ) di). 
Jo Jo 



(23) 



Thus (3 is a function of B. The solution of (23) is shown in figure 5. 
The function U(x) is now approximated for x < by 



V.(V) = C2V- 



4V£ 2 



-yV/3 v 2 ) +C 3V -*iWml, i(/Vj3 yf) 



(24) 



4V)8' 



where irj is given by (22). The magnitude of the reflection coefficient is simply the magnitude of 
the ratio C3/C2. This is determined by expanding the approximation (24) in a power series, nor- 
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FIGURE 6. Power reflection coefficient from a linearly 
stratified plasma. 




malizing it to unity at # = 0, and equating the coefficient of the x 2 term to y+(B). The result is 

(25) 



I 



Co 



y + +— - y- 






where y is given hy (20) and y_ is defined as follows: 



7-(B) = 



B 



l + -^ + 2C e — -/- + -lnj3 

4Vj8/ 2 2 2 



+ 



rVj3 



(26) 



The magnitude of the reflection coefficient thus calculated is shown in figure 6, in comparison with 
results obtained by Budden [1963] 2 by numerical integration of the differential equation. The 
agreement between the two curves is good. 



2 Private communication (unpublished numerical data). 



We are indebted to L. B. Felsen for valuable critique and stimulating discussion. The research 
which forms the basis of this paper was sponsored by the Air Force Office of Scientific Research 
of the Office of Aerospace Research, under Contract No. AF-49(638)-1402. 
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